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Abstract 

The ratio = T(rj — > 7r + 7r - 7)/r(7/ — > ir + ir~ir°) has been measured by analysing 
22 million (f> — > 777 decays collected by the KLOE experiment at DA<£NE, corre- 
sponding to an integrated luminosity of 558 pb _1 . The n — > 7r + 7r~7 proceeds both 
via the p resonant contribution, and possibly a non-resonant direct term, connected 
to the box anomaly. Our result, R^ = 0.1856 ± 0.0005 s t a t ± 0.0028 sys t, points out 
a sizable contribution of the direct term to the total width. The di-pion invariant 
mass for the r\ — > 7r + 7r~7 decay could be described in a model-independent approach 
in terms of a single free parameter, a. The determined value of the parameter a is 
a = (1.32 ± 0.08 stat l°;i° syst ±0.02 thco ) GeV- 2 . 
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1 Introduction 



The Chiral Perturbation Theory (ChPT) provides an accurate description of 
interactions and decays of light mesons [1J. The Wess-Zumino-Witten (WZW) 
term in the ChPT Lagrangian accounts for anomalous decays involving an 
odd number of pseudoscalar mesons. The triangle anomaly is responsible for 
the two-photon decays of the ir°/T)/rf mesons. Both triangle and box anoma- 
lies should contribute to the rf^ — > 7r + 7r _ 7 decays. Since the kinematic region 
of the decays is far from the chiral limit, the amplitude of the 7t + tt~ final 
state interaction has to be properly included. The decays are therefore of- 
ten described by a resonant contribution due to the p-meson exchange using 
the Vector Meson Dominance (VMD) model, and an additional contact term 
(CT), whose strength is constrained by the requirement to obtain a total con- 
tribution consistent with the WZW term in the chiral limit. In the case of 
i] — > 7r + 7r~7 the resonant p contribution is sub-dominant, making the partial 
decay width sensitive to the CT, while for the rj' — > 7r + 7r~7 decay the partial 
width is dominated by the resonance but the direct term will influence the 
shape of the di-pion invariant mass distribution. The present world average 
of the 7] — > 7r + 7r~7 partial width, T(r] — > 7r + 7r~7) = (60 ± 4) eV [2], provides 
strong evidence of the CT in the box anomaly when compared with the values 
obtained with and without the direct term, (56.3 ± 1.7) eV and (100.9 ± 
2.8) eV, respectively [3]. 

Various approaches have been used to describe the final state interaction in 
these decays: the Hidden Local Symmetry (HLS) model [3J, the chiral uni- 
tary approach [I] and the Omnes function encoding pion-pion interaction [5]. 
A model-independent approach, based on a combination of ChPT and dis- 
persion theory, has been recently proposed, where a parametrisation of the 
experimental pion vector form factor is used instead of VMD [6]. 
Recently, CLEO [7J has measured the ratio R v = T(r] — > 7r + 7r~7)/T(?7 — > 
7r + 7r~7T°) = 0.175 ± 0.007 sta t ± 0.006 syst , which differs by more than 3a from 
the average of previous measurements [HE], R v = 0.207±0.004 [TO]. We present 
a new measurement of R v with smaller statistical and systematic errors, to- 
gether with the fit of the M nn distribution according to the model-independent 
approach presented in [6]. 



2 The KLOE detector at DA$NE 



The KLOE experiment operated at the Frascati 0-factory, DA$NE, an e + e~ 
collider running at a center-of-mass energy of ~ 1020 MeV, the mass of the 
4> meson. The beams collide at a crossing angle of (tt - 0.025) rad, producing 
(f) mesons with a small momentum in the horizontal plane, = 12.5 MeV. 
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The detector consists of a large cylindrical Drift Chamber (DC), surrounded 
by a lead-scintillating fiber electromagnetic calorimeter (EMC) and a super- 
conducting coil around the EMC providing a 0.52 T field. The DC [TT], 4 m 
in diameter and 3.3 m long, has 12,582 all-stereo tungsten sense wires and 
37,746 aluminum field wires. The chamber shell is made of carbon fiber-epoxy 
composite with an internal wall of 1.1 mm thickness, the gas used is a 90% 
helium, 10% isobutane mixture. The spatial resolutions are a xy ~ 150 /im 
and o z ~ 2 mm and the momentum resolution is a(p±)/p±_ ~ 0.4%. The 
EMC [12] consists of a barrel and two endcaps, for a total of 88 modules, 
and covers 98% of the solid angle. The modules are read out at both ends 
by photomultipliers, both in amplitude and time. The readout granularity is 
~ (4.4 x 4.4) cm 2 , for a total of 2440 cells arranged in five layers. The en- 
ergy deposits are obtained from the signal amplitude, while particle position 
along fiber direction is obtained from the arrival time difference of the sig- 
nals to the photo-multipliers at the ends of each calorimeter cell. Signals of 
calorimeter cells close in time and space are grouped into clusters and the 
cluster energy E is the sum of the cell energies. The cluster time T and 
position R are energy-weighted averages. Energy and time resolutions are, 
a E /E = 5.7%/^E (GeV), and a t = 57 ps/^E (GeV) © 100 ps, respectively. 
The trigger [13J uses both calorimeter and chamber information. Data are 
then analysed by an event classification filter (EVCL), that organises data in 
different output files, according to their particle content [T4] . 



3 Event selection 



The analysis has been performed using 558 pb _1 , collected at y/s ~ 1020 MeV, 
which correspond to about 22 x 10 6 77-mesons produced. KLOE Monte Carlo 
(MC) program [T4j is used to simulate the final states produced in e + e~ colli- 
sions, taking into account machine parameters and beam-related background 
on run-by-run basis. At KLOE, the rj mesons are produced together with a 
monochromatic recoil photon, E 1 = 363 MeV, through the radiative decay 
<p — > rff. The final state under study is 7r + 7r~77 with the main background 
coming from <fi — > 7r + 7r~7r , 7r° — > 77. Another important background is 7] de- 
cay 0—7-^7—7' 7r + 7r~7r°7 — > 7r + 7r~37 with one undetected photon. In the MC 
generator the signal is simulated using a matrix element 



|M| 2 ~ fc 2 sin 2 # 



q J (M 2 - iie) 2 + M 2 P 



where, k is the photon momentum in the 77 rest frame, 9 is the angle between 
the 7r + and the photon in the di-pion rest frame, q is the momentum of both 
pions in the di-pion rest frame and V =124-(g/g ) 3 MeV with q being the 
value of q at p- meson resonance [S]. 
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After the EVCL filter, a preselection is performed, requiring at least two tracks 
with opposite charge pointing to the interaction point (IP) and at least two 
clusters in time3, not associated to any track, having energy E du > 10 MeV 
and a polar angle in the range (23° — 157°). Tracks are sorted according to 
the distance of the point of closest approach from the IP. The first two tracks 
with opposite charge are selected as pion candidates. 

3.1 77 — >- 7r + 7r~7 selection 

We require that the most energetic cluster has an energy E c i u > 250 MeV 
and we identify it as the photon (7^) recoiling against the 77 in the 0—7-777 
decay. Moreover, we ask that the 7^ is inside the calorimeter barrel (with a 
polar angle in the range 55° — 125°), to reject events with cluster split between 
barrel and endcap. To reject electrons, cuts on cluster-track association and 
identification by time of flight (TOF) are used. These cuts reject Bhabha 
scattering background and other processes with electrons in the final state. 
We exploit the — > 777 decay kinematics, to evaluate the 7^ energy: 

P M l - M l 
^ = 2 ( E ,-|p,| C0S ,) 

where, $ is the angle between 7^ and the meson momentum, p^, measured 
run by run with high accuracy using Bhabha scattering events. This allows 
us to improve the energy measurement accuracy of the recoil photon to 0.1%. 
Using and 7r-mesons momenta, we determine the direction of the photon 
(7^) from 77 decay: 

Pjr, =P<t>~ Ptt+ - Ptt- - P l4l 

the 7^ photon direction is then compared with the direction of each neutral 
cluster: Aip = ip c \ u — Lp lv (here, and in the following, the angles are evaluated 
using variables in the transverse planed). If no clusters with Aip < 8.5° are 
found, the event is rejected. The cluster with the minimum value of A(p is 
identified with 7^. In order to reject the — > tt + ?t~?t background, the angle 
between the two photons in the tt° reference frame, evaluated using the 
and the 7r-meson momenta, is calculated and it is required to be smaller than 
165°. The 7T + 7r~7 mass spectrum is shown in Fig. [lj The candidate events are 
selected requiring 539.5 MeV < M 7r+7r - 7 < 554.5 MeV. 

1 We require for each cluster \T c i u — R c i u /c\ < 5aT clu , where T c [ u is the arrival time 
at the EMC, R c i u is the distance of the cluster from the beam interaction point, 
and c is the speed of light. 

2 The azimuthal angle of the cluster is measured with an angular resolution of 6 
mrad using the position of the calorimeter cell. The polar angle is instead determined 
by the time difference of the signals at each side of the barrel and is affected by 
larger uncertainty. The use of azimuthal angle reduces the systematics. 
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Fig. 1. The tt + tt 7„ invariant mass distribution. Crosses are experimental points, 
the hashed area is the MC signal rj — > 7r + 7r~7, the filled region represents the 



total MC. Relevant background is due to 



7r + 7r 7r events and much smaller 



contribution from 0-meson decay into kaons for higher masses as well as <ft — > 177 
events for the masses below the signal peak. 



3.2 rj — > n + n n° selection 



The process — ?• 777 with rj — > 7r + 7r~7r° represents a good control sam- 
ple, having a topology similar to the signal. Moreover, in the ratio T(r] — > 
7r + 7T~j) /T(j] — > 7T + 7r~7r°) the luminosity, the <fi production cross section and 
the BR((f> — > 777) cancel out. We use the same preselection as for the rj — > 
7r + 7r~7 signal and calculate the missing four-momentum: 



P — P 

1 miss 1 q 



70 



where the variables in the formula represent the four-momenta of the <fr meson 
and of the decay products. For the 77 — > 7r + 7r~7r° sample, the missing mass 
peaks at the 7r° mass value and we select events with |M miss — M^o] < 15 MeV. 
The remaining background is rejected by an angular cut applied to the two 
photons in the tt° rest frame, ip^ > 165°. Figure [2] shows the distribution of 
the missing mass and ip^l. The two cuts select N(r] — > 7r + 7r~7r°) = 1.116 • 10 3 



events. The global selection efficiency is e 
background contamination of 0.65%. 



0.2276 ± 0.0002 with residual 
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M miss [ MeV l [degrees] 

Fig. 2. Normalisation sample 7] — > tt + tt^tt°. Left - 7r + 7r _ 70 missing mass distribution. 
Right - event distribution for the angle between prompt neutral clusters in the 7r° 
rest frame evaluated in the transverse plane, <pz?L. Crosses are experimental points, 
the hashed area is the MC r\ -4 7r + 7r~7r , the filled region represents the total MC, 
where the only relevant background contribution is due to <p — > 7r + 7r~7r events. 

4 Results 



4-1 The ratio T (77 — > tt + tt j)/T(t] — > ti + ti 7r°) 



The total selection efficiency for the signal 77 — > 7r + 7r~7 is e = 0.2131 ± 0.0004. 
In the final sample, the relative weights of signal and background components 
are evaluated with a fit to the -E m i SS — -Pmiss distribution of the 7r + 7r~7^ system, 
with the MC shapes of the remaining background and signal MC, Fig. [31 Signal 
events are counted in the range |-E m i SS — -P m iss| < 10 MeV. We find N(r] — > 
7r + 7r~7) = 204950 ± 497 events, with a background contamination at level 
of 10%. The analysis has been repeated on an independent sample selected 
without EVCL filter to evaluate any bias due to the event classification. An 
overall correction factor is used to account for data/MC difference related to 

event classification: K EVC l = iLtT'' . L + c-""° = 1-010 ± 0.009. 

^ A- — ^_4- — _0 

7T ' TT 7 7T 1 7T 7T 

Combining the results we obtain the ratio: 

= r(77 tt+tt 7) = Q lg56 ± 00Q5 ± 0.0028 

' T(?7 ->■ 71 + 71-71°) 

This result is in agreement with the recent CLEO measurement [7], while 
improving the accuracy of a factor better than three, thus confirming a smaller 
value for R v with respect to previous evaluations 09]. 

The systematic uncertainties due to analysis cuts have been evaluated by 
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Fig. 3. E m [ ss — P m iss distribution for the tt + tt 7^. The sample has been selected 
applying all the cuts described in Sec l3.ll The event counting is performed in the 
region \E miss - P miss \ < 10 MeV. 



Source of Uncertainty 


Cut Value 


Window cut 


Fractional Error on R v 


¥>77 


< 165° 


±2° 


±0.6% 




< 8.5° 


±2° 


±0.4% 


\M„+„- y -M v \ 


< 7.5 MeV 


±2 MeV 


±0.6% 


E 1 . 

mm 


> 10 MeV 


±2 MeV 


±0.1% 


CLU 


> 250 MeV 


±4 MeV 


±0.1% 


|A/ miss - M v o\ 


< 15 MeV 


±4 MeV 


±0.4% 


^77 


> 165° 


±2° 


±0.1% 


EVCL 






±0.9% 


Fit -E m iss -fmiss 






±0.6% 


Total 






1.5% 



Table 1 

Summary table of systematic uncertainties on R n . 

varying the cuts and re-evaluating the value of R v . Each cut is moved ±2<r 
with respect to the value used in the analysis, where a is the resolution on the 
reconstructed variable. The corresponding error for each source of systematics 
is reported in Table HJ The total error is taken as the quadratic sum of all of 
the contributions. 



4-2 Fit to di-pion invariant mass 



The M nn dependence of the decay amplitude has been studied in several frame- 
works. The HLS model, in particular, has been investigated in [15] and more 
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recently in [3] . In this approach, the relative strength of the CT and the reso- 
nance contribution from VMD are fixed. The model-independent approach in 
[6] , based on ChPT and dispersive analysis, does not fix this relative strength 
and parametrises the CT via a process-specific term. We use the last method to 
fit the di-pion invariant mass distribution. The function describing the partial 
width clS 8b function of = M 2 n is the following: 

= \AP{s 1T1 ,)F v {s 1 , 1T )\ r (s^) (1) 

where, A is a normalisation factor and 

r o( S nn) = 3 . 2 11 . n 3 M 3 ( M V ~ S ™) S ™ " P* 

is the simplest gauge-invariant matrix element multiplied by the phase-space 



term with j3 n = J 1 — 4M 2 /s 7r7r . F v (s n7T ) is the pion vector form factor, ap- 
proximated in the energy range of interest by the polynomial |Fy(s 7r7r )| = 
l + (2.12±0.01)s^ 7r + (2.13±0.01)4 7r + (13.80±0.14)s^, where is expressed 
in GeV 2 [6]. The P(s 7r7r ) function, a process-specific part, can be treated per- 
turbatively in the frame of ChPT, for the decay of light mesons. Taylor ex- 
pansion around s wn = gives P(s 7r7r ) = 1 + a ■ + 0(s^ w ). We fit the 
distribution by minimising the variable: 



( Ajdata y^Nbin ]\rTeo c Q \2 

x 2 = E? Un - j j j j (2) 



where, Nf ata is the content of i— th bin after background subtraction, Nj e ° 
is the content of j— th bin of the expected M nn spectrum as from Eq.(TjQ), Sj 
is the efficiency, is the smearing matrix and of = a 2 Ndata + cr^ T eo, with 
a 2 NTeo = Ef hin (Nf eo ) 2 (a^.Sf j + e 2 a 2 s ..). Figure 1 shows the measured distri- 
bution compared with results of the fit taking into account efficiency and 
smearing. 

Minimising the function in Eq. ([2]) we get 

a = (l.32 ± 0.08 stat ±°:o9 syst ± 0.02 theo ) GeV~ 2 

with x 2 /Ndf = 61/64. The theoretical error 0.02 GeV -2 accounts for uncer- 
tainty due to vector form factor parametrisation, and is determined mainly 
by the accuracy of the existing e + e~ — > tt + tt~ data. The fit is insensitive to 
the addition of a quadratic term in P(s 7r7r ). The contributions to systematic 
uncertainty on a are listed in Table [2j The value of a is in agreement with the 
result of the WASA Collaboration obtained from the fit to the 7^ spectrum 
giving a = (1.89 ± 0.25 stat ± 0.59 syst ± 0.02 thco ) GeV" 2 [16]. 
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Fig. 4. Distribution of M nn after background subtraction (black markers). His- 
togram is the fit of Eq. ([!]), corrected for acceptance and experimental resolution. 



Source of Uncertainty 


Cut Value 


Window cut 


Aq (GeV" 2 ) 


7J-+ 7f ~~ >y 


< 165° 


±2° 


+0.07/ - 0.03 


Atp 


< 8.5° 


±2° 


+0.05/ - 0.06 


\M n+7T - y -M v \ 


< 7.5 MeV 


±2 MeV 


+0.04/ - 0.04 


E" 1 . 

min 


> 10 MeV 


±2 MeV 


+0.01/ -0.04 


Total 






+0.10/ - 0.9 



Table 2 

Summary table of systematic uncertainties on a parameter. 
5 Conclusions 

Using a data sample corresponding to an integrated luminosity of 558 pb _1 , 
we select about 205000 77 — > 7T + 7T _ 7 and 1116000 rj — > 7t + 7t~7t° events from 
the — > 777 decays. We obtain the ratio of the partial widths: 

T(r] n + n--f)/T(r] ->• 71+7^71°) = 0.1856 ± 0.0005 stat ± 0.0028 syst 

in agreement with the most recent result from CLEO [7j. 
Combining our measurement with the world average value T(i] — > 7t + 7t~7t°) = 
(295 ± 16) eV [2], we find T(r] 71+77-7) = (54.7 ± 3.1) eV, which is in 
agreement with the value expected in the HLS context including the contact- 
term contribution [3]. 

We have measured the di-pion invariant mass distribution and performed a fit 
using the model-independent approach of Ref. f6j. The fit gives a = (1.32 ± 
0.08 s tat-ao9syst±0.02 tneo ) GeV 2 . 
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